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The dynamics of a thin liquid film on a solid substrate has attracted considerable interest in the past, 
allowing many applications e.g. in physics, engineering, or biomedicine (see J9) [12] HSl HS1 El E3 1 and 
the references therein). When such a film is coated with a distribution of surfactant, i.e. surface active 
agents that lower the surface tension of the liquid, the resulting gradients of surface tension induce so- 
called Marangoni stresses (e.g. E4l0 that cause the surfactant to spread. From the point of view of potential 
applications, estimates on the velocity of the spreading front and possible thinning and rupture of the film 
are among the main issues under investigation, see [ 16] and the references therein. While most part of 
the research hitherto has been dedicated to the case of insoluble surfactant (e.g., see [6, 23] and, for fur- 
ther references, see Ifl6l[l7l ). this article focuses on the actually more common case of soluble surfactants 
when there are sorptive fluxes between the surface and the bulk region immediately beneath the surface 
lTT3l [T6l [l7l [191 l23l 1241 . Sorption seems to play an important role during the early stages of the spreading. 
For instance, soluble surfactants may restrain the spreading of liquid films lfl3l or delay film rupture [23 1. 
Besides including sorptive fluxes we take into account also the effect of gravity but neglect effects of capil- 
larity and van der Waals forces. Based on the full free boundary problem for the Navier-Stokes equations 
and the mass balance equations governing the dynamics of the thin liquid and the surfactant concentrations, 
respectively, we derive a system of degenerate parabolic equations for these quantities by using a lubrication 
approximation as, e.g., in |[T2l [T3l . Cross-sectional averaging then allows one to neglect vertical compo- 
nents of the bulk surfactant concentration so that the resulting coupled system consists of one-dimensional 
equations. 

In dimensionless form the system reads: 



ABSTRACT. A thin film on a horizontal solid substrate and coated with a soluble surfactant is considered. The 
governing degenerate parabolic equations for the film height and the surfactant concentrations on the surface 
and in the bulk are derived using a lubrication approximation when gravity is taken into account. It is shown 
that the steady states are asymptotically stable. 



1. Introduction 




(1.1) 



(1.2) 



(1-3) 



for x £ (0, L) and t > subject to suitable initial conditions at t = 0: 

h(o,-) = h°, c (o, •) = c ° , r(o,-) 




(1.4) 
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and no-slip boundary conditions at x — and x = L: 

d x h = d x C = d x T = 0, x = 0, L . (1.5) 

Here, h denotes the film thickness, T denotes the surface surfactant concentration, and Cq is the (cross- 
sectionally averaged) bulk surfactant concentration. The constant f3 > indicates the degree of solubility 
so that in the limit /? — > the surfactant is highly soluble in the substrate and adsorbs weakly on the surface, 
while (5 — > oo corresponds to a surfactant accumulating preferably on the surface with weak bulk solubility. 
The surface tension er(T) is a decreasing function of the surface surfactant concentration V. The positive 
constants D and S measure molecular diffusion of surfactants on the surface and in the bulk, respectively. 
For more information on the model we refer to Section|2] 

While the derivation of models, both for insoluble and soluble surfactant, can be found in many research 
papers (the latter, however, merely in the situation when gravity is neglected), there seems to be only few 
literature investigating thin film equations with surfactant from a mathematical analytical point of view. In 
addition, to the best of our knowledge all research so far is dedicated to the insoluble surfactant case, i.e. 
when the bulk surfactant concentration Co in equations ([TT|i-([T3) (and possibly other terms) is neglected. 
For instance, in ll20l I2TI 1221 local existence results can be found for the resulting systems of degenerate 
equations for h and T. In 1 8 1 global existence of weak solutions is shown for a variant of ( |l.l| l-( fT~3l l (with 
Co neglected) including fourth order terms in h modeling capillarity effects. This fourth order degenerate 
parabolic thin film equation for h has been rigorously derived recently from Stokes flow with surface tension 
not using a lubrication approximation, cf. ifTTI . Thin films on a inclined substrate and the influence of 
various parameters on the behavior of traveling wave solutions are investigated in lfl4l [T31 l25l . We also 
refer to J3] [10] for finite volume methods and numerical simulations for thin liquid films with insoluble 
surfactant. We shall point out that research has also been dedicated to the mathematical analysis for the 
free boundary problem for two-phase flows with soluble surfactant, see fl4j[5|. Finally, we refer to the 
forthcoming paper [7| where existence of global weak solutions will be shown for the case of insoluble 
surfactant under moderate assumptions on the surface tension coefficient. 

The outline of this article is as follows. In Section[2]we provide a derivation of the governing equations 
( |l.l| >-( fT3j ). In Section [3] we begin our mathematical analysis of these equations. We establish a local 
existence result for smooth solutions and provide an energy function which shows that the only steady 
states are constants. Based on the principle of linearized stability we then show in Section |4] that steady 
states with small surfactant concentrations are locally asymptotically stable. 



2. The Physical Model 

The derivation of the governing equations for the dynamics of the liquid thin films with surfactant using 
lubrication theory can be found in many papers for various situations, see, e.g., Ifl2l [131 [T6l [T71 [191 l24l to 
name a few. The aim of this section is to provide a brief sketch of the derivation of the model equations 
( |l.l| >-( fT3j ) for a soluble surfactant including gravity that we have not found in the literature in this particular 
form. For this, we follow lfl2l[T3l to which we also refer for further details. 

We shall consider a thin film of a viscous incompressible Newtonian fluid lying on a horizontal plane. 
Initially a monolayer of soluble surfactant is deposited on the surface of the film. The ratio e of the undis- 
turbed height compared to the initial length of the film is taken to be small so that lubrication theory may 
be used. In the following, all variables are non-dimensionalized. Let x and z denote the horizontal and 
vertical coordinate, respectively. Let (u(t, x, z),w(t, x, z)) denote the corresponding velocity field of the 
fluid. The free surface of the fluid layer is at z — h(t, x) while T(t 1 x) and C(i, x, z) denote the monolayer 
concentration and the bulk concentration of the surfactant, respectively, see Figure[T] Conservation of mass 
of the fluid reads 



d x u(t, x, z) + d z w{t, x, z) = , 



(2.1) 
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FIGURE 1 . Schematic representation of a thin film with soluble surfactant. 



while the equations of momentum conservation are given by 

— d x p(t, x, z) + d z u(t, x, z) = , (2.2) 

- d z p(t, x, z) - G = , (2.3) 

where p is the pressure in the liquid and G > represents a gravitational force directed vertically down- 
wards. In ( |2.2| i, ( |2.3| ) we neglect van der Waals forces and capillarity effects, the latter assumption being 
valid if (To/5 = O(l), where S is the spreading coefficient and do is the surface tension of the uncontami- 
nated interface, see fl3l . At the bottom plane 2 = 0, no-slip boundary conditions are imposed: 

u(t, x, 0) = w(t, x, 0) = , (2.4) 

whereas at the interface z = h(t, x) the kinematic boundary condition 

dth(t, x) + u(t, x, h(t, x)jd x h(t, x) — w(t, x, h(t, x)) (2.5) 

is prescribed meaning that at the free surface the speed of the interface balances the normal component of 
the liquid velocity. Assuming zero pressure above the film and again neglecting effects of capillarity, the 
interface condition at z = h(t, x) is complemented with a normal and tangential stress condition reading 

p(t,x,h(t,x)) = (2.6) 

and 

d z u(t, x, h{t, x)) = d x a(T(t, x)) . (2.7) 

Here a(T) denotes the surface tension of the liquid related to the local surface surfactant concentration. 
Examples of constitutive relations (in dimensionless form) for the surface tension include u(T) = 1 — F or 

-3 



a(T) = (a + 1) (l - [((a + l)/a)^ 3 - l] r) 



see |6l [T2j [T3] |24) . The hydrodynamics is supplemented by mass transfer of the soluble surfactant which 
is supposed to be deposited initially at the surface of the interface as a monolayer. The dynamics of the 
surfactant surface concentration T at the free surface z — h(t, x) is governed by the advection-transport 
equation 

d t T(t, x) + d x (u(t, x, h(t, x)) T(t, x)) = D d x T(t, x) + K (C s (t, x) - T{t, x)) (2.8) 
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with D > being a non-dimensional surface diffusion coefficient (independent of the surfactant concen- 
tration). On the right-hand side of |Z8|, K (C s {t, x) - T(t, x)) with 

C s (t,x) = C(t,x,h(t,x)) (2.9) 

is the (scaled) diffusion controlled sorptive flux of surfactant between the surface and the region of the bulk 
immediately beneath the surface, where C denotes the surfactant concentration in the bulk. The constant K 
is the ratio of the time scale of the flow to the time scale of desorption. The transport equation for the bulk 
surfactant concentration reads 

d t C(t,x,z) + d x (uC)(t,x,z) + d z (wC)(t,x,z) = 8 d 2 x C{t,x, z) + 8 V 8%C(t,x, z) , (2.10) 

the positive numbers 8 and 8y in ( |2.10[ ) representing the non-dimensional horizontal and vertical bulk dif- 
fusivities independent of the surfactant concentration. Equation ( |2.10| i is supplemented with the boundary 
conditions 

8 v d z C(t,x,h{t,x)) -8d x h(t,x)d x C(t,x,h(t,x)) = -j3 K (C s {t,x) -T(t,x)) (2.11) 
at the surface z = h(t, x) and 

d z C(t,x,0)=0 (2.12) 

at the bottom z = 0. The constant (3 > appearing in ( |2.1 1[ ) is proportional to the ratio of the rate of 
adsorption to the rate of desorption, and indicates the degree of solubility. The limit f3 — > corresponds 
to the case of a surfactant with high substrate solubility only weakly adsorbing at the free surface, while in 
the limit j3 —> oo the surfactant has high surface adsorption and is only weakly soluble in the bulk. On the 
vertical boundaries x = and i = Lwe assume that there is no flux of surfactant, i.e. 

d x T{t,x) = d x C(t,x,z) = , x — 0, L . (2.13) 

Next, we shall compute u, w, and p from equations ( |2.2| >-( |2~T3"] > and then use cross-sectional averaging 
for C in order to derive the one-dimensional model equations ( |l.l[ )-( fT3] l. 
Notice first that ( |2~3| ) and ( |2~6l > yield 

p(t, x,z) — G [h(t, x) — zj . 

Together with ( |2.2| i this then implies that 

d z u(t, x, h(t, x)) — d z u(t, x,z) — G d x h(t, x) (h(t, x) — z) , 

and thus, invoking ( |2.7| i, we obtain that 

d z u(t, x, z) = d x a(T(t, x))-G d x h(t, x) (h(t, x) - z) , 

from which, in turn, we deduce by ( |2.4| > that 

u(t,x,z) = zd x a(T(t,x))-Gd x h(t,x) (zh{t,x) - y^j . (2.14) 

Integrating ( |2.1) and using ( |2.4[ ) we see that 

-h(t,x) 



d x u(t, x, z) dz = —w(t, x, h(t, x)) , (2.15) 

Jo 

whence, by ( |2.5[ ), 

/ Mt,x) \ 

d t h(t,x)+d x \ / u(t, x, z) dz = . (2.16) 



Recalling ( |2.14| > we derive the evolution equation (jTTTJ for the film thickness (for t > and < x < L): 

-h 3 d x h-- 
3 2 



d t h -d x [ ( ^h z d x h-\ h 2 d x a(T) ) = . (2.17) 
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Inserting next ( |2.14| > evaluated at z = h(t, x) into ( |2.8[ ), the transport equation for the surface surfactant 
concentration becomes (for t > and < x < L) 

d t T~d x (^h 2 Td x h~hTd x a{T)^ =Dd 2 x T + K(C s -T) . (2.18) 

We now introduce the cross-sectional average of C defined by 



1 



h(t,x) 



C (t,x):=— / C(t,x,z)dz. 
h{t,x) J 

It follows from ( |27TO] > and ( |2TT2] > that 

d t (hC Q )(t,x) = / d t C(t,x,z)dz + C(t,x,h(t,x)) d t h(t,x) 



h(t,x) 

d x (uC)(t, x, z) dz + C(t, x, h(t, x)) ( — w(t, x, h(t, x)) + d t h(t, a;)) 

h{t,x) 

d 2 C(t, x, z) dz + S v d z C(t, x, h(t, x)) . 
Owing to ( |2.5| > and ( |2.11| i, we obtain 

ph(t,x) 

d t (hC )(t,x) = d x (uC)(t 1 x,z)dz — (uC)(t,x,h(t,x)) d x h(t,x) 

Jo 

h{t,x) 

d x C(t, x, z) dz + S d x C(t, x, h(t, x)) d x h(t, x) 



- pK{C s {t,x)-T(t,x)) 

-d x ( f ( ' [(uC){t, x,z)-S d x C(t, x, z)} dz 



-/3K(C s (t,x)-T{t,x)). 

At this stage, we may assume that, within the lubrication approximation, C(t,x,z), d x C(t,x, z), and 
C s (t,x) are well approximated by Co(t, x), d x Co(t,x), and Co(t,x), respectively, and we deduce from 
( |2.14| i that, for t > and < x < L, 

*(*«,) - s, (f - i + «<w>) = -ok <c„ - r) . am 

Therefore, gathering ( |2.17| i, ( |2.18| l, and ( |2.19[ ) we obtain the model equations ( |l.l| i-( [T3| ). Notice that 

d f L , , , , „ d f L / , x 1 



, ftft, x) dx = , — / [T(t,x) + -h(t,x)C (t,x)) dx = , (2.20) 



that is, the mass of the fluid and the mass of the surfactant are conserved during time evolution, the factor 
1 //3 accounting for the degree of solubility of the surfactant in the fluid. 
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For the mathematical analysis of ( |l.l| i-( fT3] > it is convenient to introduce a new variable m := hCo//3 
replacing Co- Then the system ( |l.l| i-( fT31 l becomes (for non-vanishing h): 

d t h - d x (j h 3 d x h -^h 2 8 g a(r)j = o , 

( G 1 m \ [8 m 

d t m — d x ( — h 2 md x h h m d x a(T) + S d x m — 5 — d x h ) = —K I T 

\ 3 2 h J \ h 

d t T-d x (^h 2 rd x h-hTd x a(T) + Dd x T^j =K - T 

for x £ (0, L) and t > subject to the initial conditions at t = 0, 

h(0,-) = h°, m(0, -) = m°, T(0, •) = T° , ie(0,L) 
and no-slip boundary conditions at x = and x = L, 

d x h = d x m = d x T = , x — 0, L . 
In the following we shall focus on this version of the system ( |l.l| i-( fT3] l. 

3. Local Well-Posedness and an Energy Functional 
First we state an existence result (locally in time) for the system |Z2T]i-(|Z25|. For this we define 

O a := J ff^((0,L),M 3 )nC([0,i],(0, TO ) 3 ) , a G [0, 1] , 

where 

f {u e H 2a ((0,L),R 3 ) ; d x u = Oforz = 0,L} if a>3/4, 
H 2 N a :=H 2 N a ((Q,L),R 3 ) := { 

{ H 2a ((0,L),M. 3 ) if ae [0,3/4], 

with H 2a being the usual Sobolev spaces. Let a € (3/4, 1]. Then O a is open in Hff and a subset of 
^([(^(fLoo) 3 ). Given u := (/i, m, r), we introduce the matrices 

/i 2 



(2.21) 
(2.22) 
(2.23) 

(2.24) 
(2.25) 



a(u) := o(ft, m, T) 



/ I*. 

G h2 X 171 X 

— n m — o — o 
3 h 



a'(T) 



\ 



hmcr'(T) 



V 



G 



h 2 T D-hTa' 



:= m, T) 



/0 








KB 
IT 



KB 



\ 

-K 



K 



and we define an operator A as 

A{u)w := —d x (a{u)d x w) + b(u)w , w £ H% , u € O q . (3.1) 

Writing u° := m°, T ) equations ( |2.21| >-( [2T25] > may be recast as a quasi-linear equation in the space 
L 2 := L 2 ({0 1 L),R 3 ) of the form 

u + A(u) u = , t>0, u(0) = u° . (3.2) 
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If a is a non-increasing function, i.e. a' < 0, it readily follows that the matrix a(£i , £2 , £3) f° r G (0) 00 ) 
has only positive eigenvalues. Thus, since b(u) defines (for u G O a fixed) a bounded operator on L2, we 
infer from d Ex. 4.3.e), Thm. 4.1] that 

Ae C x -{O a ,U{H 2 N ,L 2 )) , (3.3) 

that is, A is a locally Lipschitz mapping from O a into the set of negative generators of strongly continu- 
ous analytic semigroups on L2 with domain H%, Therefore, JT] Thm. 12.1] warrants the following local 
existence result concerning equation p.2): 



Theorem 3.1. Let a € C 2 (R) be non-increasing, and let a G (3/4. 1]. Then, given any initial condition 
v? = (h , m ,r ) G Q , f/je problem \2>.2\ (and hence ( |2.21| l-( p.25| l) admits a unique (strictly) positive 
strong solution 

u = (h, m, r) G C([0, T), o a ) n C7 1 ((0, T), L 2 ) n c((o, T),H%) 

with maximal time of existence T G (0, 00]. 

Notice that, thanks to the positivity of h, Co := f3m/h is well-defined having the same regularity prop- 
erties as m and (h, Co, T) solves ( |l.l| i-( fT31 l. Notice further that the positivity of solutions is built into the 
set O a . 

Next we shall derive an energy functional for the system ( |2.21[ )-( |2~25| l. 

Proposition 3.2. Consider a non-increasing function a G C 2 (K), a G (3/4, 1], and an initial condition 
vP = (h°, m°, r ) G Q . 77ze corresponding solution u — (h, m, T) to ( |2.21| l-( p.25| l given by Theorem 
satisfies 



3.7 



<9, 



-\£h \d x a(r)\ 2 dx-£(^ d x h - ^ ft 1 /* a^r) j dx 



/W 2 



da; 

2 



n 



for t G [0, T), the function <j> being such that 

(t>"(r) r = -o-'(r) > , r > . (3.4) 



Observe that each term on the right-hand side is non-positive, in particular due to the monotonicity j3A\ 
ofcf>'. 
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Proof. We deduce from p3l) -< |ZB) that, for t € [0, T), 



d_ 
dt 



cj>{T) + -h<t> 



j3m 
~h~ 



cj> l {T)d t Y + 4> 1 



dx 

/3 m 

IT 



d t m + - 



P to 



/3 m , ( P m 



h 



h 



d f h > dx 



J <f>"(T)d x r (-^h 2 Td x h + hTd x a(T)-Dd x rj dx 

' /3m 



A" 



K 



1 



P m 
~h~ 

(3 m 
~h~ 



-T da 



(3m 
~h~ 

j3 m 
~h~ 



G , n r m 

/T m + £ — 

3 /i 



1 



- r dx 



Recall that (13. 4b ensures 



/3ra^ /3m 

~7T 



9 r 



/3 m 
~h~ 



h 3 d x h 

3 



d x h — (5 c^to + -mh d x a(T) > dx 



1 



hf d x a{T) } dx 



b"(T)Td x T = -d x a(T) and hd x 



(3m 
~h 



= (3d x m 



j3m 



d x h . 



Hence, noticing that the last integral of the right-hand side of the above equality cancels with the first and 
the last term of the third integral, we have 



da; 



[ L h\d x a(T)\ 2 dx-D [ L c/>"(r)\d x T\ 2 dx-^ f 
Jo Jo P Jo 



j3 m 

~h~ 



9y 



f3 m 



dx 



K 



4>'{T) - 



j3 m 



/3m 



dx+ ~ [ h 2 d x hd x a(T)dx 
2 Jo 



From equation ( |2.21[ ) it follows that 



* G [\ 2 dx = - G2 
dt 2 / n 3 



I h 3 \d x h\ 2 dx+ - / h 2 d x hd x a(T) dx 
Jo 2 J 



which, added to the previous equality and after completing the square, yields the stated energy equality. □ 



An interesting consequence of Proposition 3.2 is that it provides a complete description of the (smooth 
positive) stationary solutions to ( 2.2 1 [ >-( 2T25] > when a is decreasing. Indeed, observe that, if u is a (smooth 
positive) stationary solution to ( |2.21| i-( |2~25] l, the left-hand side of the energy equality vanishes and so does 
each term of the right-hand side since they are all nonnegative. It is then straightforward to establish the 
following result: 

Corollary 3.3. Suppose that a € C 2 (]R) is strictly decreasing. Then the only (smooth positive) steady 
states to ( |2.21[ )-( |2.25[ ) are of the form (ft*, to*, r*) with constants ft* > and m* — ft*T*//3 > 0. 

The next section is dedicated to the local asymptotic stability of the steady states (ft*, to*, T*) for small 
values of the surfactant concentrations, i.e. for small values of to* and r*. Let us mention at this point that 
some constraints have to be satisfied by the initial data because of the conservation of the mass of the fluid 
and the mass of the surfactant ( |2.20| i. 
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4. Asymptotic Stability 



To prove stability of steady states we use the same notation as in the previous section and write equations 
( |25T) -( |Z25) in the form Q with u = (h, m, T) and u° = (ft , m°, r°). 

Let ft* > be a fixed constant and, given a sufficiently small number 77* > 0, set 

ft* 



Notice that /3to* = ft*r* and thus, u» := (ft*,TO*,r*) is an equilibrium of ( |3.2| i. Motivated by ( |2.20[ ) we 
then introduce a projection P e £(£2) n C(H^), i.e. P 2 = P, by setting 



Pw := 



ft - (ft) 



ft* +/3 





V ft* +/3 



where 



(m + T) 



to + r) 



1 



= (h,m,T) € L 2 , 



denotes the mean value of a given function / 6 L 2 ((0, £)). Then both spaces L 2 and P^ decompose into 



L 2 = PL 2 8 (1 - P)L 2 , 4 = PP^ © (1 - P)P 



Also note that 



(l-P)A(v)w = 0, (1-P)b(v)w = 0, we P^, weOi 



that is, A(v)w £ PL 2 and b(v)w £ PL 2 for any v £ 0\ and w G P^r. Thus, if u is any strong solution to 



([3T2J on [0, 00) with initial condition u° = (ft , to , T°) <E H% satisfying (ft ) = ft* and (to + T° 
then necessarily 



(4.1) 
1 1< 

V* 



(1 - P)u(t) = (1 - P)u° 



t > 



Therefore, we have 



u{t) = v{t) + u* e PHfj 8 (1 - P)H 2 N , t>0, 



and v := Pu solves 

v + A*v — (A* — A(it* + u)|pjj 2 ) v ~ H u * + v ) u * —'■ F{ v ) ) ^ > 



u(0) = P?i° = u° - u* =: u° 



(4.2) 
(4.3) 



with 



A*w := A(u^)\ PH ^ w + B*\ PH 2 w , w £ PH% . 
Here, the operator P* , given by 

/ \ 

P/3to* 



P* w = 



hi 
K /9 m* 



, w = (h, m, r) e p 1 , 



(4.4) 



is the Frechet derivative of [z H> b(z)u*] : Oi/ 2 L 2 at z = u*. 

Lemma 4.1. 77ie operator A*, considered as an unbounded operator in PL 2 with domain PJJ^, belongs 
to r hL{PH%, PL 2 ), that is, —A* is the generator of a strongly continuous analytic semigroup on PL 2 . 



II) 
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Proof. Noticing that A(u* ) |(i_p)jj 2 = ar, d recalling ( |3.3| > and ( |4. we may interpret A(u* ) as a matrix 
operator 

' A (u*)\phI ° X 











G n(PHff © (1 - P)H 2 N ,PL 2 © (1 - P)L 2 



We then conclude from E I.Cor. 1.6.3] that A(u*)| PH 2 G n(PH%,PL 2 ). NowB.I^ G C{H 2 Nl PH 2 N ) 
can be considered as a perturbation, see [2, 1.Thm. 1.3.1]. □ 

Next we show that — A* has negative type for small values of 77*. 

Lemma 4.2. There are numbers e :— e(h*) > and loq := u>o(h*) > such that the spectrum er(— A*) of 
—A* is contained in the half-plane [Re A < — cjo] provided that < 77* < e. 

Proof. Let w° = (ft , m°, T°) G be arbitrary and let iu(t) := e -4 " 4 * w°, rj > 0, be the unique solution 
in PL2 to the linear equation 

w + A*w = 0, t>0, w(0) = w°. 

Writing w = (ft, to, T) we have 

(h(t)) = , (m(t) + r(t)> = , t > , (4.5) 
and, by definition of A* and (|4.4|>, 



9tm — 9 X ( — ft* to* G^ft — - ft.* to* <r'(r*) c? x r + 5 d x m — 5 — — d x h 



ihh-n, [~hld x h-~hlo-'(r*)d x T] =0 



G 



1 



- — m- KT - K (5 — h = {) , 
ft* hi 



G 



d t T ~d x [^ hi r* 9 x /i + (£> - ft, r* a'(T,)) d x T 



I ^A m + KT + Kp r ^h = 0, 
h* hi 



for x G (0, L) and t > subject to 

d x h = d x m = d x Y = , x = 0,L, t>0 
Given q > we define the symmetric matrix 

/ 



b q (h*,m*,T*) := 



q 4ki 

3 

/3G, 2 £/3to* 

— — ft* to* — — 

6 2 ft* 



'G 2 5/3 to* G 3 g 2 , \ 



TO* ft* cr'(r*) 



G 



V^^r* - jhla'{r,) -^TO*ft*d'(r*) Dft* -ft2r*d'(r*) y 



/3 



Then, multiplying the equations satisfied by ft, to, and T by qh, fim, and ft*T, respectively, we derive the 
equality 

1 d 

K 

ftl 



(?||fc||2 + 0IMl2 + Mr||!) + (b q {K,m.,T t )d x w 



(),. 



K fim 

+ — \\l3m- KTg + — ^ (ft|ft* r - pm) 2 = 



(4.6) 
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where || • || 2 and (-|-)2 denote the norm and inner product in L 2 - Next, since 



we have 

det(6 9 (fc*,0,0) - A) = {6 p- A) 







o 





\hla'{Q)\ 





Dh* 



A - 



DL A 



!GD^-^«<7'(0) 2 



so that all eigenvalues of b q (h*, 0, 0) are positive provided that 

16GD 

° <q< 3^W 

Thus, for q satisfying this condition, there is some e := e(/i*) > such that the matrix b q (h*, to*, T*) is 
positive definite for < 77* < e by continuity and the definitions of to* and T* . From ( |4.6| l we then derive 
that 



1 d 
2di 



(g|Wl2 + j8|M|l + Mr||l)+/*(||^| 



\d x m\\ 2 2 + \\d x T\\l) 



K 

/?.* 



(4.7) 



for some number fj, := /i(/i*) > 0. Recalling ( |4.5| l we may apply Poincare's inequality to deduce that 

1 

2 

Hence, writing 

1 , „ , , „ 1 



il^mlll 



\d x T\\l> -\\d x m + d x T\\l>c\\m + T\\l and \\d x h\\j > c\\h\\ 2 



(4.8) 



to = - — — (/3 TO 
/3 + K v 



/i» r + h, (to + r)) 



Y = j^(p{m + T)-{pm-K T)) 



we obtain from ( |4.8| > that 

IMI2 + lirill < c (||/3m - Jur||| + ||a s m||5 + Rru 2 ) . 

Plugging this into ( |4.7| > and making e (and hence m* < 77* < £ on the right-hand side of ( |4.7| i) smaller if 
necessary, we deduce that 

\^(q Ml + P \\m\\l + K ||r|||) + Mo (9 + /3 ||m||l + ^ ||r|||) < 

with a sufficiently small number (j,q := fio(h>*) > 0- This, in turn, readily implies that w = (h, to, T) has 
exponential decay, i.e., 



-tA,„..0 



W\\ n < Me- Uot \\vP\U , t>0, 



for some numbers M := M(7i*) > 1 and := wqC 1 *) > 0. This yields the assertion. 



□ 



Finally, note that the function F on the right-hand side of ( |4.2[ ) is defined on any sufficiently small 
neighborhood V of zero in PH^ and that, since b(u* )u* = 0, 

F G C 2 ~ (V, PL 2 ) , F(0) = , F'(0) = , (4.9) 

with F' denoting the Frechet derivative of F. The smallness condition on V is needed to guarantee that the 
first component of 7i* + v does not vanish so that 6(71* + v) is well-defined and we may in particular choose 



the zero neighborhood V in PH^ such that 



< min{/i*, to*, T*} , veV 



(4.10) 
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Lemma 4. 1 Lemma |4~2] and ( |4.9[ > are then the ingredients to apply the principle of linearized stability [18 
Thm. 9.1.2] to ( |4.2| >, ( |4.3j ). We thus conclude that given any cj G (0, cjo) there exists R > such that for 
any initial value v° belonging to a ball in PH]^ centered at zero of sufficiently small radius, the problem 
( p~2] >, ( |43| ) admits a unique global solution v € C 1 ([0, oo), PL 2 ) H C([0, oo), Pf/^) with e V for 
i > and 

Mt)\\pH% + \Wt)\\pL 2 < Re- Ut \\v°\\ PH% , t > . 

Therefore, u := [t H + it*] is the unique solution to ( |3.2) , Notice that ( |4.10| i and u(i) G V for t > 
warrant the positivity of u. 

Summarizing we have shown the following result on the local asymptotic stability of steady states for 
( |2.21| i-( |2~25] l with small surfactant concentration: 

Theorem 4.3. Let a G C 2 (M.) be decreasing. Let h* > be arbitrary. Then there exist numbers 
e := e(h*) > 0, u> := uj(h^) > 0, and M := M(/i*) > such that for 

with rj^ 6 (0, e) and any initial value (h , m° , T°) G with (h ) = /i* and (m + T ) = 77* satisfying 
the smallness condition 

\\h° - h4 H 2 +\\m°- m4 H 2 + ||r° -r*||fla < e , 

there is a unique global positive solution (h, m, T) G C 1 ([0, 00), L%) n C([0, 00), H^) to p.21| i-( |2~25] l. 
This solution satisfies 

\\h(t)-h,\\ H 2 + \\m{t)-m,\\ H 2 + ||r(<)-r,|| H 2 

< Me~ ut (\\h° ~ h4 H 2 + \\m° - m4 H 2 +\\T° - T^) 

for t > 0. 
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